P-wave Pairing in Two-Component Fermi System with Unequal Population near 

Feshbach Resonance 
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We explore p-wave pairing in a single-channel two- component Fermi system with unequal popu- 
lation near Feshbach resonance. Our analytical and numerical study reveal a rich superfluid (SF) 
ground state structure as a function of imbalance. In addition to the state A±i oc Yi±i, a multitude 
of "mixed" SF states formed of linear combinations of Yim's give global energy minimum under a 
phase stability condition; these states exhibit variation in energy with the relative phase between 
the constituent gap amplitudes. States with local energy minimum are also obtained. We provide 
a geometric representation of the states. A T=0 polarization vs. p-wave coupling phase diagram 
is constructed across the BEC-BCS regimes. With increased polarization, the global minimum SF 
state may undergo a quantum phase transition to the local minimum SF state. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.50.-s 
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Introduction: Discovery of s-wave paired fermion con- 
densates in cold atoms [l| subjected to s-wave Feshbach 
resonance (FR) have led to explorations in several fas- 
cinating directions. One is the study of superfluidity in 
asymmetrical Fermi systems, a subject of intense exper- 
imental and theoretical research [3, l3|, Ij]. Another is 
the prospect for attaining fermion condensates with pairs 
having non-zero relative orbital angular momentum. Re- 
cent observations of p-wave FR [5|, [y, |7|, |8| in ^Li and 
^'^K have raised the possibility of observing p-wave su- 
perfluidity in cold atoms, leading to renewed theoretical 
interest. P-wave superfluidity in the BEC-BCS crossover 
region has been studied using fermion-boson models ap- 
plied to single-component Fermi gas (where FR occurs in 
the same hyperfine state) [9|, llO| , and using fermion-only 
models at T=0 both in two-component Fermi gas (where 
FR occurs between different hyperfine states) [11[ and in 
single-component case [12|. Ohashi [l3| considered both 
cases at finite-T. While there have been fervent theoreti- 
cal efforts [31, |j| on pairing in asymmetrical Fermi systems 
subject to s-wave FR, work on the corresponding p-wave 
case is lacking. 

In this paper, within a fermion-only model, we study p- 
wave superfluidity near a p-wave FR in a two-component 
Fermi system with unequal populations in two hyperfine 
states. It is conceivable that in future, population imbal- 
anced systems under p-wave FR would be realized similar 
to the ones subject to s-wave FR. We hope that this work 
will motivate experimental efforts. 

It has been pointed out [7|, [llj that unlike liquid ^He, 
pairing interaction in cold atoms is highly anisotropic in 
"spin"-space. ("spin" referring to hyperfine states). For 
example, in ^Li, when the hyperfine pair |l/2, —1/2 > is 
at resonance, the pairs |l/2, 1/2 > and | — 1/2, —1/2 > 
would not be. So, in studying pairing in this two- 
component system, intra-spin interactions need not be 
considered. Pairs in p-wave superfluids with unlike 
"spin" components can however have different orbital an- 



gular momentum content, so the gap parameters A/m ex 
i ~ 1 spherical harmonics Yim {m = ±1, 0). 

Our analytic and numerical considerations of the pop- 
ulation imbalanced system lead to several interesting re- 
sults and predictions: (i) We find a rich superfluid ground 
state (GS) structure: The p-wave state Ai (hence A_i 
by symmetry) give a GS global minimum energy. Inter- 
estingly, under a specific condition involving the relative 
phases of the three pairing amplitudes A^'s, a multi- 
tude of "mixed" states of the form oAq -I- 6Ai -I- cA_i 
(Aim = Am here) are degenerate with A-ti. In addition, 
we find states with local energy minimum, (ii) We pro- 
vide a geometric representation of the p-wave superfluid 
states (Fig 1): The states exhibiting global minimum lie 
on a "semicircle" formed by the intersection of the surface 
of the sphere formed by Ai , A_i , Aq with a plane defined 
by |Ai|-|-|A_i| = const, (iii) We obtain polarization (P) 
vs p-wave coupling phase diagram at T=0. (Fig 2). The 
superfluid phase SFl comprises of states with global min- 
imum, while SF2 of states with the local minimum. We 
also find a region of phase separation PS between SF2 
and the normal phase N. In this two-component system, 
PS persists onto full polarization, (iv) These raise the in- 
triguing possibility for a T^O quantum phase transition 
from SFl to SF2 at finite polarization (Fig 2). Addi- 
tionally, transitions at finite-T would also occur, (v) In 
the limit P — > 0, we find that the ground state structure 
of P ^_0j_ is preserved, hence richer than that obtained 
earlier 
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Model: We consider a two-component Fermi system 
with unequal "spin" (t,i) population, but equal masses 
for the unlike fermions. We take the interaction between 
unlike fermions to be isotropic in orbital space. Following 
the rationale above, interactions between like fermions 
are taken to be zero. Since we adjust self-consistently the 
chemical potential with the strength and sign of the cou- 
pling, in our fermion-only model, molecules would nat- 
urally appear as 2-fermion bound states. The pairing 



Hamiltonian is then given by: 

kcr 

+ X! ^kk'Ck+q/2|C_k+q/2i'^-k'+q/2iCk'+q/2T (1) 
kk'q 

where Cko- (ck^.) is the annihilation (creation) operator for 
a fermion with momentum k, kinetic energy ^ko-^ek— Mo-, 
and spin a (=Tj i); /^<t is the chemical potential of each 
component, and eit^—fi^k^ /2m. 

We consider condensate pairs with zero center-of-mass 
momentum, (q=0). TC is mean- field (MF) decoupled via 
the "spin" -triplet (5* = l,ms = 0) pairing gap function 
A||(k) = A(k)=-X;k' Vkk'(c-k'TCk'i) giving: 



kcr 



[A(k)4^cLkx+H.c. 
-^|A(k)|V^kk 



E 

k 



(2) 



To obtain the variational ground state energy of the MF 
Hamiltonian ^ we use the equation-of-motion method 
[j| for the imaginary-time (r = it) normal {G^a' ) and 
anomalous {Fa-a-i ) Green's functions: 

drGaa'iKT) = -S{T)Saa'-£.kaGaa'{KT) 

+ A_,,(k)F_,,,(k,r), (3) 

drF„„,{k,T) = ^^kaFaa'(k,T) 

+ A;_,(k)G_..,(k,r). (4) 

Eqs. ^ and ([4]) are Fourier transformed with 
T— >iLj„=i^ and dr^'—v, where iuJn~{2n-\-\)'KT are 
Matsubara frequencies. This gives Gaa^k^v) = 
-(C-k-.+^)/[(a.-i^)(e-k-.+i')+|A,._,(k)|2] 

and i^,_,(k,i.)-A;^,(k)/[(ek-. - i^)(e-k. + 1^) + 
|A,_,(k)|2]=._F_„,(k,i.). 

We separate the radial and angular parts of the inter- 
action potential: Vkk' = (47r/3)E™^fcfc'>"i,m(k)n!m(k') 
(k=(0, 0)), and take for V^.^.' a separable form conve- 
nient for analytic calculations: Vf.f.i = Xw{k)w{k'). 
Using (c_k'j,Ck'|) = F|*|(k',T=0^) and Fourier 
transforming, we obtain for the gap function, 

A(k)=-(l/A:BT)^^k,Vkk'-FxV^»e"°^- ^ith the 
separable form of Vkk' , the gap function can be written 
as: A(k) = Y.m^mw{k)Y^M = w{k)A{k), where 
A™ = -il/kBT)j:^^,Xwik')Y,Uk^)F^^{k',,.)e'^^^\ 
With i^||(fc,i/) known, we finally obtain the equations 
for the gap amplitudes (m ~ 0, ±1), 

Am = -Xj2w{k)Y*^{k)A{k)[np{E^)--np{E+)]/{2Ek) 

k 

(5) 



where E^ = -h±Ek,h= (^h - ^k^)/2 = {p^- - ^ J/2, 
i?fe = [e^ + |A2(k)|]i/2, with ^fc=(efcT+ai)/2=efc-Ai and 
/i=(/i|+((i|)/2; np{Ef. ) are Fermi functions. We can now 
average ^ and calculate the variational ground-state en- 
ergy Eg- At T=0 and for /i>0 this is given by: 



Y,{{Ek-h)e{-Ek+h)}--y.\^ 
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(6) 



where the regularized interaction g is given by 
l/.g=l/A+(l/(27rft)3) / w^{k)(fk/2ek. 

Ground State Analysis: First, we carry out an 
analytic study of the p-wave superfluid ground state 
structure. We assume |A(k)| <C Ek, and expand 
Ek in powers of |A(k)| keeping terms to 4th order: 
Ek ~ |6|[H-|A(k)|V2e2-|A(k)|V8ea. Substituting 
this expression for Ek in ([6]) and converting momentum 
sums to integrals we obtain a simple representa- 
tion for Eg = aJ\ACk)\^dn + f3j\ACk)\^dn + -f, 
where a = {l/{27rhf) J^^^f^w^{k)Pdk/2ek + 
il/{2nhf)J^^^^wHk)k^dk{l/2ek-imk) - l/.g, 
(3 = {l/{27:hf)J^^^^wHk)k^dk/8\^k\' and 

7 = (47r/(27r/i)^) /g ^^(^fc-/i)fc^dfc. Existence of a 
superfluid phase requires a<0, otherwise minimization 
of Eg forces the gap to vanishes. The polarization 
(P = (?i| — ni)/{ni + 7i|)) dependence of Eg is con- 
tained in a,/3,7, which depend on h, fj,,'w{k), while a 
alone depends explicitly on the coupling g. Owing to 
these dependences the gap parameters are sensitive to 
changes in P. 

For fixed ^ and /i, and Ek given above, the ground state 
energy Eg is obtained by minimizing ([6]) with respect to 
the pairing amplitudes A^. Thus we arrive at an analytic 
expression for the ground state energy [14| : 



Eg = -^+1 + 2P (^|Ao|^ + \A,f + |A_i|2 + ^ j 

-H/3(|Aop-2|Ai||A_i|)2+4/3(l-i)|AonAi||A_i|(7) 

where t=cos0, with 9 = cos^^ (AoA^f, A^A-i) = 
(200— '/'I— 0-i); ^m's are the phases associated with the 
gap amplitudes A^. From (O, we find that for a stable 
p-wave superfluid phase to exist, the following conditions 
have to be satisfied simultaneously 



(a) |Ao|2+|Ai|2+|A_i|2 = -a/4/3 
(6) |Ao|2-2|Ai||A_i| =0 
(c) (l-t)|AonAi||A_i| =0 



(8) 



giving the GS global minimum energy 

Sg = -aV8/3 + 7- (9) 

Conditions (a) and (b) together represent a semicircle 




FIG. 1: (Color online) Geometric representation of p- 
wave states in two-component population imbalanced system. 
States (e.g. A, B, C) exhibiting global energy minimum lie 
on the semicircle formed by the intersection of sphere surface 
with plane (see text). For clarity, plane containing semicir- 
cle is also shown separately (green). O (north pole) and M 
(on Ai-A_i plane) are states with local energy minimum. 
Am(»Ti = 0,±l)'s have been normalized to the sphere radius 
v/-a/4/3; a < 0. 



formed by the intersection of the surface of a sphere 
of radius R, R^ = —a/4/3, with a plane defined by 
All + |A_i|=_R. All points on the semicircle must nec- 
essarily satisfy the phase constraint (c). They repre- 
sent a multitude of "mixed" superfluid states of the form 
aAo -I- &Ai + cA_i . A 3D geometric representation of the 
states is shown in Fig. 1. The GS global energy minimum 
condition is satisfied only for i = 1 in Eq. 8c, i.e. when 
the relative phase 9 = ((^o— '/'i) — (0-i— 0o)) = ^nn. 
For the states A±i, (points A and B on the semicircle 
in Fig 1) this is always satisfied. For the other states 
on this curve, containing varying admixtures of A^'s, 
Eg fluctuates with 0, and attains the global minimum 
only for 9 = 2mT (Fig. 3a); the fluctuation ampli- 
tude depends on the specific state chosen; e.g. C at 
the top of the projected figure has the maximum am- 
plitude variation, but still lies lower than the local min- 
imum we discuss below. t=l corresponds to a paral- 
lel orientation of vectors AqA^ and AqA_i; this gives 
the same global minimum i?G=— a^/8/3 + 7 for three 
particular cases: (A) |Ao| = |A_i|=0 and |Aip=-a/4/3; 
(B) |Ao| = |Ai|=0and|A_ip=-a/4/3; (C) \Ao\^=-a/8f3 
and |Ai 1^ = 1 A_ip= — a/16/3. These are shown as points 
A, B,C respectively in Fig. 1. In addition to these, we also 
find the existence of local minimum with £'q=— a^/12/3 + 



and |Aip=|A_i|^=— a/12/3; these are shown as O and 
M respectively in Fig. 1. Values for the GS energy are 
completely determined by a, /3, and 7. In evaluating 
these, we take the N-SR [l5| form: w(fc)=fcofc/(fco + k^), 
where ko is a cut-off momentum. 

In the limit of zero polarization, P ^ 0, the GS en- 
ergy coefficient 7 — > 0, while a < 0,/3 > remain finite. 
Thus the GS structure of the polarized case is preserved 
for zero polarization, with the same expressions for the 
energy global and local minima (different numerical val- 
ues). The states A and B (equivalent to A by symmetry) 
above and in Fig. 1 correspond to the finding of Ref. [llj . 
Our work reveals an additional set of "mixed" states on 
the intersecting semicircle (Fig. 1) that also have the 
same global minimum for specific values of the relative 
phase discussed above. 

Numerical Calculations: Guided by our analysis of the 
p-wave superfluid ground state, we carry out a detailed 
numerical study. For fixed number densities, we solve 
sclf-consistcntly three gap equations (Eq. 5) and two 
number equations given by: 



J2'\''ka 



{4aCka) 



k 



G^,{k,T^O-), 



(10) 



7 for: (i) |Aor=-a/4/3 and |Ai 



=0; (ii) lAohO 



This gives us the gap amplitudes A^, and the chemical 
potentials fia-- Using these we obtain the ground state 
energy Eq from Eq. 6, as well as the coefficients a, (3, 7 
that determine Eg given by Eq. 7. The agreement be- 
tween our numerical Eq and our analytical Eg (based 
on an expansion of Ek) is good to within 15%. 

We construct a polarization (P) - coupling (l/kpOt) 
phase diagram in BEC-BCS crossover regime (Fig. 2), 
enforcing the stability conditions: that the stability ma- 
trix {d^ Eg /d Amid Amj) is positive definite; and that 
Sp/Sh > 0. In Fig. 2, SFl denotes a stable superfluid 
phase, corresponding to the states on the 'semicircle' in 
Fig. 1 that produce GS global minimum, i.e. with rel- 
ative phase 9=2m: among the gap parameters. With 
increased polarization, at T=0, SFl becomes unstable, 
and gives way to the superfluid phase SF2, correspond- 
ing to states with the local minimum discussed above. In 
Fig. 2, states Ai and Aq were chosen in SFl and SF2 re- 
spectively; other choices of allowed global and local min- 
imum states give qualitatively the same phase diagram. 
At even larger polarizations, SF2 becomes unstable to 
phase separation, PS. SFl, SF2, PS occupy a relatively 
much narrower part of the phase diagram on the BCS side 
compared to the BEG side. In our two-component system 
with inter-species interaction, PS persists into full polar- 
ization, P=l. This is reasonable because at P=l, the 
system is essentially a one-component system in which 
the absence of minority species atoms makes inter-species 
interaction inoperative. Such a system can exhibit super- 
fluidity only under the effect of intra-species interactions. 
Our results suggest several interesting possibilities: A 
r=0 quantum phase transition from SFl to SF2 driven 
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FIG. 2: (Color online) Polarization P vs p-wave coupling 
— l/(A;|.at) phase diagram of two-component one-channel 
Fermi system with p-wave pairing. Shown are normal, N, and 
p-wave superfluid phases, SFl and SF2, as well phase separa- 
tion, PS. Unitarity limit is shown by the dashed vertical line. 
Inset: Calculated chemical potential vs coupling across BCS 
and BEC regimes for P = 0.1 (solid line); 0.5 (dashed line). 



by polarization may be expected. Since the energies of 
SF2 and the continuum of "mixed" states for 6 ^ 2mT 
lie higher than the GS global minimum, they may be ac- 
cessed at finite-T. Also, each p-wave state on the semicir- 
cle (Fig. 1) (generally formed of a linear combination of 
the A,„'s except for the "pure" A±i states A,B), is char- 
acterized by a relative phase 2mT < < {2n+2)TT. Hence, 
for each there is a multitude of 6'-depcndent states (Fig. 
3a) with varying Eq that may be interesting to probe 
with experiments sensitive to such relative phases. The 
inset in Fig. 2 shows the calculated behavior of chemical 
potential /i across the BEC-BCS regime. It deviates sig- 
nificantly from the Fermi energy in a wider region around 
FR even on the BCS side, and drops much more rapidly 
to negative values on the BEC side compared to the s- 
wave FR case. For sufficiently weak coupling in the BCS 
regime, /i approaches Fermi energy ep. 

Figs. 3b, 3c show the variation of Eq with coupling 
g {=25kpat/8TT) at fixed polarizations P, and with P 
for fixed g. Eq is normalized to the Fermi energy ep 
= h'^kp/2m, with 2fc|, = k^,-^ + kp^. For a given P, 
Eg becomes less negative as g approaches unitarity; the 
trend is more noticeable for smaller P's. For a given 
g, Eg lie higher for smaller P, presumably due to the 
lower majority-species band becoming progressively more 
occupied with increasing P, thereby lowering Eg with 
increasing polarization. For small P, Eg becomes less 
negative with increasing g. This trend is reversed for 
P> 0.3. 

Conclusions: We believe we have presented new re- 



FIG. 3: (Color online) (a) Calculated ground state en- 
ergy (scaled to global minimum) vs. relative phase an- 
gle 9 (see text). Solid curve with maximal fluctuations 
of Eg corresponds to state C in Fig 1, (|Aop=— a/8/9, 
\Aif = \A-if^~a/mi3). The dash and dot-dash curves 
correspond to two other states on the arc in Fig. 1: 
|Aop=-a/11.3/?, |Aip=-a/76.3/3, |A_ip=-a/6.73/3, and 
|Ao|^=-a/16/?, |Ai|^ = -a/186/3, |A-i|^=-a/5.49/3 respec- 
tively, (b) Calculated ground state energy vs. coupling g for 
different polarizations, P. (c) Calculated ground state energy 
vs. polarization P for different coupling strengths g. 



suits and made several predictions for two-component 
population imbalanced Fermi systems across the entire 
BCS-BEC regimes. Future experiments in cold fcrmions 
near p-wave FR may be able to provide tests of some of 
these predictions. Insight into the nature of the orbital 
part of our superfluid states may be gained from measur- 
ing the angular dependence of momentum distributions; 
from molecular spectroscopy using light radiation; or pos- 
sibly measurements of zero sound attenuation. The work 
could also be of interest to other two-component Fermi 
systems where p-wave intra-species couplings are small 
or negligible. 
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